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Abstract. The linear stability of accretion disks is revisited. The governing equations are expanded asymptotically 
and solved to first order in the expansion parameter e defined by the ratio of the disk's vertical thickness to its 
radial extent. Algebraically growing solutions are found for global perturbations on the radial accretion flow of 
thin inviscid compressible Keplerian disks. The algebraic temporal behavior is exhibited in the vertical velocities 
and the thermodynamic variables and has the form t sin Qot locally in the disk where Qo is the Keplerian rotation 
rate. The physical implications and relations to the Solberg-Hoiland stability criteria are discussed. 
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1. Introduction 

It is well known that the pure Keplerian shear flow (KSF 
in this letter) is linearly stable, i.e. no exponentially grow- 
ing solutions are known. This stability analysis is imple- 
mented in accretion disks (AD) to also imply their lin- 
ear stability. However, the governing equations for AD 
entail the existence of meridional circulation in the disk 
(Kluzniak & Kita, 2000, Regev & Gitelman, 2002) and 
here we are interested in how weak meridional flow affects 
the linear stability of ADs. 

A classical approximation in ADs is to integrate 
over the thin vertical extent of the disk. However, since 
the dominant energy transfer is in this direction, this 
may result in a loss of some salient physical features. 
Furthermore, it is also known that the meridional circula- 
tion has a significant component in this direction. Hence, 
the present work includes the vertical dimension explicitly. 

To the point: the possibility that transient growth 
(TG) phenomena may play a significant role in astro- 
physical disks has gained recent attention with the work 
of Iaonnou and Kakouris (2001) and Yecko (2004). The 
essence of the TG phenomenon is that significant growth 
of an initial perturbation develops (up to orders of mag- 
nitude) before final decay, (cf. for a review see Schmid 
& Henningson, 2001). Indeed, linear stability analysis of 
KSFs show that asymptotic stability is promoted because 
(i) the Rossby number of the flow is ~ 1 and, (ii) there 
are no inflection points in the KSF profile (Balbus, 2003), 
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except in the star-disk boundary layer (Bertout & Regev, 
1996). 

Even though linear disturbances decay in the long run, 
the fate and character of short time transient responses 
(and its nonlinear consequences) in an AD may be signif- 
icant and deserves deeper exploration. Two observations 
provide circumstantial support to the idea of TG of KSFs : 
(a) the analogy to laboratory flows like Couette- Taylor ex- 
periments (Richard & Zahn, 1999, Longaretti, 2002) which 
have many common features with KSFs and (b) linear 
stability analysis of shearing flows in localized sections 
of KSFs show significant TG (Chagelelishvili et al, 2003, 
Tevzadze et al, 2003, Yecko, 2004, Mukhopadhyay et al., 
2004, Afshordi et al., 2004) leading to significant nonlinear 
ramifications (e.g. 2D simulations of Umurhan & Regev, 
2004). Also, Iaonnou & Kakouris (2001) demonstrate TG 
for global disk perturbations but their considerations were 
restricted to two-dimensional (radial-azimuthal) distur- 
bances of an incompressible KSF. Studies of viscous com- 
pressible boundary layer shear flows have been shown by 
Hanifi & Henningson (1998) to have solutions which grow 
like t in the inviscid limit. What is most fascinating is 
that these algebraically growing solutions are present in 
the limit where the streamwise perturbations are absent. 
These trends beget the question: is there some sort of anal- 
ogous compressible TG effect in accretion disks in the in- 
viscid limit? 

In this work we show that algebraically growing so- 
lutions indeed do exist for global perturbations of a geo- 
metrically thin AD. The steady state is assumed to be a 
general barotrope and we consider relatively fast adiabatic 
perturbations. The equations studied here emerge from an 
asymptotic analysis of the equations appropriate for such 
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a disk environment (Regev, 1983, Kluzniak & Kita, 2000, 
Regev & Gitclman, 2002) in the zero viscous limit. 

2. Equations and general solution 

We assume that what governs the evolution of the geomet- 
rically thin ADs are the Euler equations (inviscid flow). By 
'geometrically thin' we mean that 



e = H/R < 1, 



where the vertical and radial scales are H and R re- 
spectively. The governing equations are asymptotically ex- 
panded in e and this is detailed in Appendix A. We as- 
sume the expansion converges fast for sufficiently small 
e. The physical advantage of the asymptotic expansion is 
that the procedure draws out of the governing equations 
terms with equal physical importance, namely the terms 
with the same e power have the same order of magnitude 
effect on the flow. Further, we assume here that the per- 
turbations can be treated adiabatically, namely we restrict 
the discussion to time scales shorter than the dissipation 
timcscales in the AD. The three momentum conservation 
equations governing the evolution of disturbances in cylin- 
drical coordinates are (to the lowest non vanishing power 
in e), 



d t u[ 
rd t n' 2 
Podtv' 2 



-u'^d r r 2 n Ql 



-d z P 2 - p 2 g(r,z), 



(2) 
(3) 
(4) 



where the dynamical radial (r), azimuthal (</>) and vertical 
(z) velocities are denoted by u^fl^v^ respectively. The 
mass continuity and energy conservation equations are 



dtp' 2 = -^d r paru[ - d z p v' 2 , 

d t P' 2 = -u' 1 d r P -v' 2 d z P ~ 1 P { 1 -d r 



d z v' 2 ) . 



(5) 
(6) 



The above are nondimensionalized according to the proce- 
dure described in Appendix A. We stress that the undis- 
turbed state contains the possibility of steady meridional 
circulation only for non-zero viscosities. The pressure and 
density disturbances are P 2 ,p' 2 while their steady state 
configurations are Po, pg. The vertical component of grav- 
ity (which points towards z = 0) is expressed generally 
as a function of r and z; for thin rotationally supported 
Keplerian disks it is given by, g{z,r) = z/r 3 , up to lowest 
nontrivial order in the e expansion. The steady rotation 
rate for Keplerian flow is given as Qq = r~ 3 / 2 . The ratio 
of specific heats is 7. 

In the spirit of this general discussion, it is assumed 
that the steady state thermodynamic quantities behave 
as barotropes which means that there is a unique function 
relating P = Po(po). For general barotropes the spatially 
dependent barotropic index, n(r,z), is given by = 
j j p ° . The steady state relationship satisfies the verti- 
cal hydrostatic equilibrium, namely, d z P$ = —g(z,r)po. 
It should be noted that integration of this equation in- 
troduces an arbitrary function of the variable r which is 



usually taken to be the height of the disk, h(r), which is 
assumed to be well behaved. Finally, we define the adia- 
batic sound speed as c 2 = 'fPo/po- Once the form of the 
barotrope is known then so is the sound speed. Eqs.J2EJ) 
may be reduced to the simple PDE, 

(d 2 + n 2 e ) u[ = 0, j = r 2 n , nl = ^d rJ 2 . (7) 

If as an initial condition we have an arbitrary meridional 
flow field and zero azimuthal flow perturbation, i.e. at t = 



(1) 0, u[ = u(r,z), n' 2 



u(r, z) cos fl e t, 



0, then the solution to (J7J is 
Q,' 2 =(ismti e t, (8) 



with fl = — u(r, z)/2r 1 / 2 . The epicyclic frequency Q e is 
equal to flo in KSF's. Note that because f2 e is a function 
of r, inspection of the solution at this order reveals its 
inseperability. Combining (] ll(il) reveals the single PDE, 

{d 2 + C) p v' 2 = d z (u'^rPo + 2SigWi) + f cW«i,(9) 
where the differential operator C is given as 



[ft 



[9 



n-\~l — 717 
n + 1 



)] 



~^gd z - c 2 d 2 z 



(10) 



We notice immediately that the evolution of v' 2 is explicitly 
driven by the solution obtained for u[ . With p' 2 = P 2 = 
at t = we have the solutions, 



M^h sin tt e t + A4 0) cosfi e t, 
p {1) t cos n e t + p {0) smQ e t, 
P {1 Hcosri P t + P (0) sinQ.t. 



Pov 2 

P' 2 = P (1) icosft e < + P (0) sinft e i. (11) 
The structure vector M = (M^\ M^ ] ) satisfies 
[£-n2]AfW=#i i = 0,l (12) 
where the source vector 3>, the driver, is given by: 

* = {-fa,-(4> a + 2n e Mi 1] )}, (13) 

in which 

0o = d z {ud r P Q ) + d z (-yP ~d r ru) + ^d r p ru, 

4>i = - (d z -fP u + gp u)d r D, e . (14) 

The remaining structure functions are given to be 



P 
p<°> 
p(i) 

n P p^ 



1 



.p(l) 



,(1) 



r ru - 



£ (d z M^ - Po ud r a 

\d r rpm + 9,Mf + p 
gMW - jPod^M^/po + jPoudrfl 
ud r P Q - 7P0 (7^ 



For d r Q e 7^ the expression for <f>i is in general not zero. 
The above family of solutions is completely determined 
once the initial perturbation in the radial (meridional) ve- 
locity, u (or equivalently rh = pott), is given (in contrast to 
a density or a pressure perturbation). Next, the structure 
vector M must be determined subject to the boundary 
conditions. The classical assumption is the vanishing of 
the Lagrangian pressure perturbation on the moving sur- 
face at z s = ±h(r,t) (e.g. Korycansky & Pringle, 1995). 
We also study the homogeneous solutions of (J5J. These 
acoustic modes (i) oscillate stably and (ii) have no fre- 
quencies which are resonant with sinf2 e i (for n constant). 
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3. Nature of the solution 

Classical stability analysis of linear systems assumes solu- 
tions of the simple exp(iwi) form and searches for a disper- 
sion relation for ui. Here we were able to find a complete 
solution (i.e. ~ t$mVl e t) that cannot be represented by 
a finite number of exponential functions. The basic rea- 
son for this behavior lies in the structure of the properties 
of the e expansion as reflected in the the set I|7l9f) and 
revealed in 112fl . The equation for the first component is 
linear in the perturbation. The second equation is also 
linear however, the second equation has as a source the 
solution to the first equation. This is a system of stag- 
gered linear equations where at each level the solutions to 
the previous level appear as source terms. 

The physical effect responsible for the algebraic growth 
can be traced from the driver vector <!> to arise from the 
basic shear, i.e. <9 r £! e . Solutions growing with t depend on 
<pi. This growth vanishes if there is no shear. 

The main physical consequence of the solution is an 
algebraic growth on the shortest time scale in the problem, 
namely the local dynamical time scale r = 1 / rZ e . The 
growth continues at least till the next order terms in e 
become important. The exact point of saturation (in this 
inviscid limit) is under investigation now. However it is 
clear that the smaller e = H/R is the longer the growth 
will continue. The effect is inherent to the thinness of the 
disk (or equivalently how cold it is). Also, the solution 
found here is in the limit of vanishing viscosity. In reality, 
viscous terms would enter at some stage to cause a decay. 

The existence of algebraically growing solutions is in- 
sensitive to the global entropy gradient, i.e. whether or 
not it is stable to buoyancy oscillations, or to the sign of 
the Raylcigh criterion, i.e. d r j. This may then give the im- 
pression that the solution found here violates the Solberg- 
Hoiland criteria (SHc) (e.g. Tassoul, 2000, Rudiger et al., 
2002) which predicts the instability of linear infinitesimal 
axisymmetric disturbances of steady rotating gas flows. 
The criteria pertain to the dynamics of infintesimal per- 
turbations of the steady solution of the full Euler equa- 
tions IjA.ll) and the adiabatic condition (|A.2|I without any 
assumption about the size of e. We refer to this pertur- 
bation procedure as classical-linearization (CL). The re- 
sulting equations are subsequently analyzed and, following 
some simplifying assumptions, imply the SHc. 

On the other hand, the approach we have taken here is 
different in philosophy: we instead seek to develop a finite 
amplitude non-linear solution to the problem of dynam- 
ical disturbances of the full equations. In this approach 
we start with the apriori assumption that e is small. This 
is followed by assuming that solutions of the equations, 
both steady and dynamical, may be e expanded. Solutions 
are determined by iteratively solving the resulting equa- 
tions at each order of e. Unlike the CL procedure, in which 
disturbances arc infinitesimal, disturbances here are intro- 
duced with amplitudes comparable with the characteristic 
parameter of the system, namely e. It is in this sense that 
these perturbations are considered finite-amplitude. 



Because the two approaches yield different evolution 
operators (see below) it is not surprising that the algebraic 
growth we discover is not predicted by the SHc. The SHc 
are not valid for finite amplitude perturbations as these 
but, instead, they are valid for infinitesimal ones. 

To put this in somewhat more concrete terms let 
Vl = (u L ,il L ,v L , p L , P L ) T , represent a column vector 
of the perturbed quantities resulting from the CL proce- 
dure. In CL the governing equations for the evolution of 
the perturbed quantities are, 

d t v L + vv L = o, p=(d l 6 )' (15) 

with the following operator submatrices defined: 

/ o - 2 2\ ( ft £M 

H= i, ' , L= d zPs 0, 

Kl * ]r J \d z p s + ~fP s d z o o J 

in which j r — d r j, and, 

\ d r P s + jP s ~o r r / v 7 

with g r — ——d r P s . The quantities P s and p s are the 
steady state pressure and density solutions to (|A. 1|) for 
arbitrary e. The SHc refers to the perturbation matrix V 
of the linearized system. No reference to the e parameter is 
made. We can go one step further and expand the matrix 
operator V into powers of e and when doing so we discover 
that the O e operator, since it is proportional to e 2 , can 
be neglected for thin accretion disks (see below). 

By contrast, having started from (A.l) and expanded 
in e, we arrived upon a different system: one which has 
taken advantage of the extreme geometry of the sys- 
tem (i.e. e <C 1) but also one which has bypassed the 
assumption of linearization and all of its host implica- 
tions. In particular, the finite-amplitude expansion pro- 
cedure implemented has to lowest order lead to lll'Kil) . 
which is expressed compactly, dtV' + VqV' = 0, with 
V = (%, fl 2 , v 2 , p 2 , P2) 7 '■ This equation is similar to l|15|l 
in which Y 1 — > V and where Vq is V. The differences 
between these are (i) P s — * Pq and po — > p s and, (ii) O e 
is replaced by O, a 2 x 3 zero matrix. Evidentally V and 
Vq differ in quality due to the absence of O e - 

The algebraic growth solution discovered results from 
this finite-amplitude approach. The obvious question to 
ask is: how do the two approximations relate to each 
other? The approximation carried out here is not re- 
stricted to infinitesimal perturbations but instead is one 
that exploits the geometrical nature of the object under 
discussion. It is therefore clear that the implications of the 
SHc analysis and the analysis presented in this work refer 
to two disjoint perturbation spaces which evolve according 
to different operators. 

It is easily seen that the work performed by the solu- 
tion JSJ) per cycle (27r/f2 e ) is equal to 2 / K~(PQu{d r Sl e /VL 2 e ). 
Evidentally this work is zero if the shear vanishes. It is 
mainly this extra work which drives the algebraic growth. 
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4. Physical meaning and consequences 



Appendix A: Equations, scalings and expansion 



ADs arc subject to continuous noise. The consequences of 
the effect discovered here is a fast rise in the vertical di- 
rection followed, presumably, by decay. The process takes 
place continuously. The disk is never 'quiet' . The abil- 
ity to observe such a rise depends on the radial extent of 
the initial perturbation, ft is clear that the perturbation 
should be observed to propagate from small to large radii. 

The main result of this work shows that both the ver- 
tical velocities and the density/pressure fluctuations show 
algebraic growth proportional to t sin n e t where f2 e is the 
disk rotation rate while the radial velocity (and m) show 
no growth in this approximation - again a result of the 
staggered nature of the system of asymptotic expansion. 
Physically the result holds so long as t < 1/e for otherwise 
the asymptotic orderings loses their validity. Finally the 
results apply to any rotation law provided the object has 
e < 1. 
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We consider the equations of circumstellar flow under the 
influence of a central gravitating point source and subject 
axisymmetric perturbations. Self gravity is ignored. The 
scalings here are identical to those used in the germinal 
investigations of the these problems, (e.g. Regev, 1983, 
Kluzniak & Kita, 2000, and Regev & Gitelman, 2002). 
The Keplerian rotation speed (fifc) is set by the Keplerian 
value at the fiducial radius, i.e. tt 2 ~ GM/R? with corre- 
sponding Keplerian speed V% ~ GMj R. The gas pressure 
is given by P = c 2 p where c is the speed of sound and the 
pressure scale is given by P ~ pc 2 . 

Since H ~ c 2 /g, it also happens that e 2 = c 2 /V 2 . The 
choice for the temporal scale is the sound crossing time in 
the vertical direction, which happens to be equivalent to 
the rotation time at the given radius R. We assume that 
the radial (u) and vertical (v) velocities are scaled by the 
sound speeds and the rotation rate (f2) is scaled by flk- 
Substituting all of the above into the governing equations 
yields the following, 



ed t u - 
pd t fl 



e 2 ud r u 



ev d z u - fTr = -e z ?-d r P 



i 

77 



1 



e^d r r 2 n + pvd z n = } 



dtv + eud r v + vd z v 



-\d z P 



77 



dtp + t\d r rpu + d z pv = 0. 



(A 



We assume here short time scale dynamic perturbations so 
that it is safe to assume adiabatic perturbations, dS/dt = 
0, where S is the specific entropy of the gas given by 
CylnP/p 1 , where 7 is the usual ratio of specific heats 
and C v is the specific heat at constant volume. Under 
these assumptions the energy equation becomes: 



(d t + eud r + vd z ) P + jP (e±cU + d z v) = 0. 



(A.2) 



We consider the following asymptotic expansions of the 
solutions 

n = O (r) + e 2 [n a (r, z) + Q,' 2 (r, «,*)] + ••■ 

u = e [ui (r, z) + u[ (r, z, t)} -\ 

v = e 2 [v 2 (r,z) + v' 2 (r,z,t)] ^ 

P = P (r, z) + e 2 [P 2 (r, z) + P 2 '(r, *,*)] + ••■ 
p = p (r,z) + e 2 [p 2 (r,z) + p , 2 (r,z,t)} + --- (A.3) 

To lowest order we find the standard thin disk solutions 
which do not depend at all on any viscous mechanism: 

n Q = r- 3 / 2 , d z P = -p g(r,z), g(r,z) = z/r 3 , (A.4) 

together with barotropic equation of state. The equa- 
tions governing the dynamics of the temporally evolv- 
ing quantities (i.e. p 2 , P 2 , v' 2 , u[, Q! 2 ) are quoted in the 
body of the text. The solutions to the steady quantities 
P2, Pi, v 2 , Mi, ft 2 are not sought here but are detailed in 
Umurhan et al (2005, in preparation). We do note that 
in the absence of viscosity Vi and Ui vanish to all orders 
unless we externally impose some amount of radial mass 
flux into the system. 



